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Condensation of atom pairs with finite total momentum is expected in a portion of the phase 
diagram of a two-component fermionic cold-atom system. This unusual condensate can be identified 
by detecting the exotic higher Landau level(HLL) vortex lattice states it can form when rotated. 
With this motivation, we have solved the linearized gap equations of a polarized cold atom system 
in a Landau level basis to predict experimental circumstances under which HLL vortex lattice states 
occur. 
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I. INTRODUCTION 

Polarized two-component fermion systems tend toward 
finite pair momentum condensates because of the Fermi 
radius mismatch between majority and minority compo- 
nents. In superconductors, electron spin-polarization can 
be induced by the application of an external field or by 
proximity coupling to a ferromagnet. Finite-momentum 
Cooper pair condensates in spin-polarized superconduc- 
tors, Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states, 
were first proposed in the early 1960's PJQ- One im- 
portant consequence of finite-momentum pairing in an 
isolated superconductor is a spatially inhomogeneous or- 
der parameter. There have been many efforts in various 
solid state systems to detect this exotic state, including 
recent ones 0-EI- but its definitive identification has re- 
mained elusive. The disorder that is inevitably present 
in a solid state system may have played a role in the ab- 
sence of a conclusive FFLO state identification in studies 
of spin-polarized superconductors. 

Experimental progress U El @ with fermionic cold- 
atom systems has given rise to a new strategy for re- 
alizing the FFLO state or the related Sarma state 
and has stimulated a great deal of theoretical activity 

PIlllIipiIllIlplIipiIll^l^l^l^P- 

The tunability of the interaction between atoms via a 

Feshbach resonance [51 EH has made it possible to in- 
crease the strength of fermion pairing and has even made 
the BEC-BCS crossover [22J |28|, |2i| experimentally ac- 
cessible. On the Bose-Einstein condensate (BEC) side 
of a Feshbach resonance fermionic atoms form bosonic 
molecules which condense at low temperatures. On the 
Bardeen-Cooper-Schrieffer (BCS) side, the effective at- 
tractive interaction between fermion atoms leads to BCS- 
typ e pairing. In between lies the so-called unitarity limit 
|30j in which no weakly-interacting particle description 
applies. 

Easy control over the population of two hypcrfinc 
states in a trapped atom cloud, makes cold-atom systems 
a promising candidate for FFLO state realization. The 
FFLO state competes EILlEEQEOllllE 
El El El El El with a number of other states, in- 



cluding in cold atom systems states with phase separated 
regions that are respectively unpolarized and unpaired. 
The FFLO state is expected to occur on the BCS side of 
the BEC-BCS crossover, at temperatures and pressures 
close to the normal/superfluid phase boundary. Popula- 
tion imbalance in cold atoms plays essentially the same 
role as a Zeeman or exchange held in a superconduc- 
tors since pairing is dependent on energy measured from 
the Fermi energy for each species of fermion. In both 
cases the Fermi radius of the majority species exceeds 
the Fermi radius of the minority species and pairs at the 
Fermi energy necessarily have non-zero total momentum. 

One of the most obvious signatures of superfluidity in 
fermionic cold-atom systems is the appearance of vortices 
and vortex lattices when the system is rotated |3l|]. In- 
deed recent experiments |j| have observed vortex-lattice 
structures in fermionic cold-atom systems close to the 
BEC-BCS crossover region. For this reason an obvious 
potential signature of an FFLO state is the appearance 
of the exotic vortex-lattice structures they are expected 
to form miHHI]. FFL0 vortex lattices can be wildly 
different from the usual hexagonal Abrikosov vortex lat- 
tice. The structure of the vortex lattice is determined 
mainly |3l l33l l34| by the Landau level index of its con- 
densed fermion pairs; the Abrikosov lattice forms when 
the Landau level index j = 0, which is the closest ap- 
proximation to zero-total-momentum pairing allowed in a 
system that has come to equilibrium in a rotating frame. 
FFLO states in the absence of rotation can imply j > 
Fermion pair condensation in rotated systems. Vortices 
have been observed in systems with population imbal- 
ance @, but so far no unusual vortex structures have 
been observed. (This could be due to the fact that these 
experiments realize the gapless Sarma phase |'2 and an- 
other reason could be that the FFLO state is predicted 
by weak-coupling theory while all experiments are in the 
unitary limit.) 

With this motivation, we report on a study of 
the polarization and interaction strength regime over 
which non-zero j pairing is expected in a rotating two- 
component Fermion system. We consider only the BCS 
side of the Feshbach resonance, on which FFLO physics 
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FIG. 1: Low energy pairings for population balanced and 
unbalanced systems. Shaded regions indicate participating 
states for the low energy parings in fc-space. Q is the total 
momentum of the pairs, which is for balanced systems and 
equal to the difference between Fermi wavevectors in unbal- 
anced systems. 



occurs. We consider three-dimensional systems for the 
sake of definiteness, although two-dimensional systems 
could also be interesting experimentally. Working in the 
co-rotating reference frame, rotation is equivalent to an 
external magnetic field and a reduction in radial con- 
finement strength. All our explicit calculations are for 
a uniform three-dimensional system and do not account 
for confinement. In typical experiments the atomic Lan- 
dau level splitting, equal to 2MI where VL is the rotation 
frequency, is much smaller than the Fermi energy. In 
this limit the Landau level index of the condensate could 
be determined by finding the optimal pairing wavevector 
on the BCS superfluid/normal phase boundary in the 
absence of rotation and using semiclassical quantization 
to add rotation to the condensate effective action. Here 
we use a fully quantum-mechanical approach, including 
Landau quantization even at the level of the underlying 
unpaired fermions. This approach is still relatively easy, 
partly because of the short-range of the atom-atom at- 
tractive effective interaction, and has of the advantage 
of determining the condensate Landau level index more 
accurately, and allows us to comment on the rapid rota- 
tion regime which might be approached experimentally 
in the future. From now on we use the language of the co- 
rotating frame so that the atoms experience an effective 
field with cyclotron frequency Q c = 2f2. 

Pairing is most effective when the states to be paired 
are as close to the Fermi energy as possible. When there 
is no population imbalance, pairs formed from electrons 
with opposite momentum (zero total momentum) are 
abundant at low energies as illustrated schematically in 
Fig. ^ For unbalanced populations the lowest energy 
pairs have total momentum equal to the difference be- 
tween Fermi wavevectors. In systems with an orbital 
magnetic field linear momentum is not a good quan- 
tum number, but the motion of a pair can still be sep- 
arated into center-of-mass and relative motion degrees- 
of-freedom. In a magnetic field, momentum space col- 
lapses into Landau levels whose degeneracy is illustrated 
in Fig. [5] by partitioning of momentum space into equal 
area segments centered on hrt c (N + 1/2). A pair of elec- 
trons with given Landau level indices N and N' has fi- 
nite quantum amplitudes for all center of mass Landau 




FIG. 2: Degeneracy of Landau levels. States between two 
dotted circles collapse into the solid circle. All the areas be- 
tween two adjacent dotted circles are the same and solid cir- 
cles have radii given by h 2 k 2 /2m — hQ c (N +1/2). The arrows 
show the maximum and minimum momentum differences be- 
tween particles in LL N — 1 and N = 2, which correspond 
qualitatively to the maximum and minimum of the COM mo- 
mentum. 



level indices from to N + N' which correspond closely 
to the distribution of center of mass (COM) kinetic en- 
ergy values that would be obtained by averaging over the 
corresponding regions of momentum space illustrated in 
Fig. |2 These quantum probability amplitudes are the 
key ingredient in the linearized gap equations discussed 
below. We derive linearized gap equations which implic- 
itly define the critical temperature for a phase transition 
from the normal to the superfluid state for each COM LL 
and determine the phase boundaries in parameter space. 
If excited COM LL's have a higher critical temperature 
than the lowest-lying COM LL, this signals the occur- 
rence of exotic vortex lattice states and of FFLO states 
in the unrotated system. In Fig. [3] the phase diagram 
is shown for a non-rotating homogeneous system. The 
maximum value of the exchange field (or difference be- 
tween normal state chemical potentials) for which pairing 
still occurs is given approximately by H = Ao/y/2, where 
Ao is the BCS gap parameter at zero exchange field and 
zero temperature. Beyond this so-called Clogston limit 
[3(| the BCS state is no longer stable. The FFLO state 
is expected to occur in this region of the phase diagram. 

The rest of the paper is organized as follows. In Sec. ITU 
we derive COM Landau level index dependent linearized 
gap equations for the critical temperature of the rotating 
system. The numerical solution of this equation is pre- 
sented in Sec. IIIII We finish in Sec. IIVI with a discussion 
of our results, and present our conclusions. We postpone 
to this section a discussion of the competition between 
phase separated states and FFLO states, which is an is- 
sue for cold atoms but not for electrons in a solid because 
of long-range repulsive Coulomb interactions. 



3 



Normal 



0.707 
0.704 




N=50, M=50 



N=40, M=60 



FIG. 3: BCS theory phase diagram for FFLO and BCS states 
as calculated, for example, in Ref.[3j|. Here, H is the ratio 
of the Zeeman energy (or normal state chemical potential dif- 
ference) to the zero-field energy gap. The dotted line marks 
the Clogston limit where the energies of normal and the zero- 
pairing momentum BCS state are identical. The FFLO state 
occurs near the boundary between normal and BCS states. 
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II. LINEARIZED GAP EQUATIONS 

In this section we derive the linearized gap equation for 
condensation of Fcrmion pairs with a definite COM Lan- 
dau Level (LL) index. We first consider the transforma- 
tion between individual particle and COM and relative 
states for two rotating atoms and then use this to derive 
the gap equations, which are implicit equations for the 
critical temperatures of each COM LL index channel. 



A. Unitary Transformation 

To consider the pairing instability of a normal Fermi 
gas, we first turn our attention to the description of scat- 
tering between two atoms in a rotating reference frame. 
The rotation is represented by considering the atoms to 
be particles with unit charge in an effective homogeneous 
orbital magnetic field. The Hamiltonian for two particles 
is 

h = ^ H*V P1 - A(n)) 2 + (-ihV r2 - A(r 2 )) 2 
4(- ffl V R - 2 A(R)) 2 + l(- !ft V r -^ 



(1) 

where M — 2m, fi — m/2, R = (ri + r 2 )/2 and r = 
ri — r 2 . The vector potential A(r) is defined by V x 
A(r) = 2mQ, z where 51 is the angular rotation frequency 
of the system and we assume that the rotation is around 
the z-axis. In the Landau gauge, A(r) = (0, 2mflx, 0) 
and the individual atom eigenfunctions with eigenvalues 
hn(2N + 1) are given by 



where <^>Ar(r) is the one-dimensional harmonic oscillator 
eigenfunction and the effective magnetic length Ib is de- 
fined by h 2 /ml B — 2hQ. The eigenfunctions are labeled 
by the momenta in y and z directions, and by the LL 
index N. The eigenfunctions for the COM and relative 
coordinates are the same, except that the effective mag- 
netic lengths are now Ir — lg/y/2 and l r — \21b- In 
terms of ladder operators, 

h = Tiil c (a\ai + a\a 2 + 1) = Ti£l c (a} R a,R + a\a r + 1) , (3) 

where cti = (Ib/ v2h) (7r,> - ini >y ), Tti = iTiV t - A(r l ), 
o-b = ifli + CI2) /v2) a r — ifli — CL%)/y2, and Ml c = 
H 2 /ml B = 2MI. The ladder operators can then be used 
to derive [37| an explicit expression for the unitary trans- 
formation between individual particle and COM and rel- 
ative two-atom states: 

(ri,r 2 |JV, fa >y , ki, z ;M, k 2 , y , k 2 , z ) 

N+M 

= E Bf M (R,r\j,K y ,K z ;N + M-j,k y ,k z ) , (4) 



where 



and 



K y = ki >y + k 2 , y , K z = k ljZ + k 2 . z , 

k v = ( k i-,v - k 2.,y)/2 , k z = (ki, g - k 2 , z )/2 , 



B 



NM 



j\(N + M — j)\N\M\ 



nl/2 



2N+M 



{-) 



M—m 



(j - m)\(N + m - j)\(M - m)\m\ 



(5) 



4>n{Xi 



■hjDKLyLz) 1 ' 2 , (2) 



It follows that Bj is the probability amplitude for two 
atoms in LLs N and M, respectively to have COM LL j 
and the relative motion LL N + M — j. When TV = M, 
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N',k' y ,k[,^,v' M'-k' r -kl,i,a-v' 

FIG. 5: Ladder diagrams to be summed for scattering func- 
tion r. 

\Bj fM \ 2 has maxima for j = and j = N+M. However, 
if N ^ M, \B^ M \ 2 can have a maximum for intermedi- 
ate j, which means that for two atoms in different LLs, 
the most probable COM LL can be different from zero 
or N + M as shown in Fig. 0] The smooth envelope ap- 
parent in these figures is simply the zero-field probability 
distribution of the COM kinetic energies given the Fermi 
momenta of two individual particles. The COM energy is 
maximum for parallel momentum and minimum for op- 
positely oriented individual particle momenta. This coef- 
ficient plays an important role in determining the pairing 
COM LL in condensed states. 



B. Bethe-Salpeter Equations 

The pairing instability in a Fermi gas is signaled b 
a divergence of the many-body scattering function |38j 



which we approximate using the Bethe-Salpeter equa- 
tion summarized by the finite-temperature Feynman di- 
agrams illustrated in Fig. We consider a system con- 
sisting of two hyperfine species denoted by ] and J, . For 
definitcncss we assume that the two species have the 
same energy spectrum but allow for different densities 
and therefore different chemical potentials. Population 
imbalance is relatively easy to achieve experimentally and 
the life-time of each hyperfine state is long enough com- 
pared to experimental time scales to justify the use of 
equilibrium statistical mechanics with separate particle 
reservoirs for the two species. The many-body scatter- 
ing function is calculated by summing the ladder dia- 
grams ^^^(see Fig. El . Generalizing the calculations 
of Ref.[ 37] to three dimensions from two and we find that 
the total two-particle scattering function can be written 
as a sum over different COM Landau level index chan- 
nels: 



T(N, M, k y ,k z ;N', M', k' y , k' x ;iu) 

= Bf M Bf' M ' lj {N, M, k y , k z ,N', M', k' y , k' z ;iu)(6) 

3 



where the partial scattering function for COM LL j 



(N + M- j, k y , k z \V\N' + M' - j, k' y , k' z ) 
E E \Bf M "\ 2 (N + M - j, k v ,k z \V\N" + M" - j, k'^, K 

N",M" k%,k'> 

x K N ,,,M'',k'j{iu>)^ j (N' , ,M'',k''k'^N',M',k' k' z .;icj) . 



(7) 



In Eq. 



1 - f(^N,k^) ~ f{€M,-k*,0 



1 



£N,k z = hfl c [ N+- 



h 2 kl 



2m 



(8) 
(9) 
(10) 



and /(£) is the Fermi distribution function. In the case 
of a delta-function interaction V(r) = — Vo8(r) we have 
that 



(N + M- j, k y ,k z \V\N' + M' - j, k' y X) 

~Vo4 l N+M-j(kylr)4'N<+M'-j(kylr){l/LyL z ) 



(11) 



where <jf N is the one-dimensional harmonic oscillator 
eigenfunction in relative coordinates. Using this property 



and the orthogonality of the relative motion harmonic os- 
cillator wavefunctions we find that 



7j-(JV, M, k y ,k z ;N',M',k' y ,k' z ;iu;) 



To 



LyL z 



x 1 



Vo 



B N",M",k'J 



E K N n M „ k ,,(iuj) 



B 



n" m" 



(12) 
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C. T c Equation 

As mentioned before, the instability of the normal state 
due to pairing is signaled by the divergence of the many 
body scattering function T{iu> = 0), and therefore a di- 
verging -fj(iuj — 0) means that pairs with COM LL j are 
unstable to condensation. This instability condition for 
the scattering function is equivalent to the lineariz ed g ap 
equation which defines the critical temperature (38j in 
mean-field theory. (In the mean-field-theory for the or- 
dered state |40j the order parameter can be expressed 
in terms of partial contributions from each COM LL 
channel. When the order parameter is small the vari- 
ous channels decouple and the partial contribution from 
a given channel vanishes at the same point at which the 
normal state partial scattering function diverges.) From 
Eq. I|12|) , we get an implicit equation for the critical tem- 
perature Tl for each COM LL j, which reads 

1 _ 1 1 - f(€N,k z ,l) ~ f(tjM,-k z ,j) \qNM\ 2 

(13) 



Unlike the BCS superconductors, for which retarded 
phonon-mediated attractive interactions have a natural 
ultraviolet cut-off, there is no cut-off in this equation and 
the summation is over all states. Hence, as it stands, 
this equation diverges, because of the assumption of a 
(5-function interaction. To remove this divergence, we 
need to recognize that the true atom-atom interaction 
is short-ranged compared to relevant atomic wavelengths 
but not a ^-function. Using the exact relation between 
scattering length and interaction strength (see Eq. I|A.5|) 
in the appendix) we remove the interaction strength Vq 
by renormalizing to the scattering length (3^ in the 
equation and obtain convergent sums over intermediate 
states. The equation for T-! then becomes 



k F oCl s 



2nk F Q 2— 1 

e u N,M ' 



dk. 



D 



NM I 



£N,k z + £M,-k z 



LB 



NM I 



(14) 



where kFo is the Fermi wavevector of the unpolarized 
system without rotation. The left-hand side of Eq. (|14|> 
is experimentally measurable. We determine T c as a 
function of l/fci?o fl sc by solving this implicit equation 
combined with implicit equations for the temperature- 
dependent chemical potentials fi a 



(15) 



N,k y ,k z 



where n a is the density of atoms in hyperfine state a, 
and V is the total volume of the system. In the next 
section we present numerical results obtained by solving 
these equations. 



III. NUMERICAL RESULTS 

We calculate T-! for each COM LL j for various rota- 
tion frequencies, interaction strengths and polarizations 
and in this way determine the phase boundaries in the 
parameter space spanned by hfl c , a sc and the polariza- 
tion. We fix the total density of the system n to t and 
used the polarization p as a parameter. The polarization 
is defined by 



P : 



71 i 



(16) 



where n-\ is the density of the majority species and rij_ is 
the density of the minority species. Hence, the density 
of atoms in species a [a = +1 (—1) corresponds to f (j)] 
is given by 



I + up 



■ "tot 



(17) 



The relationship between T-! and interaction strength 
is illustrated in Fig. The true critical temperature for 
the system is the largest value of . 



T c = max {T^ } 



(18) 



At weak rotation [Fig. (a)] , the transition temperature 
T c for zero polarization shows the usual behavior |4l) 
T c (x exp(— l/kFoa S c) and the highest T-! is for the j = 
channel regardless of the interaction strength. In this 
circumstance we expect the system will have a standard 
Abrikosov vortex lattice. The critical temperature de- 
creases as polarization increases and superfluidity is sup- 
pressed above some critical polarization. It is more easily 
suppressed at weak interaction. FFLO states, which cor- 
respond to nonzero j, occur at strong interaction and 
high polarization. We emphasize that these states will 
have very distinct vortex lattices, more open than 
the hexagonal Abrikosov lattices and qualitatively dif- 
ferent for each value of j. It should be quite obvious 
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-1 / k vr a 

FO sc 



(a) Weak rotation 




-1 / k m a 

FO sc 

(b) Medium rotation 



phase boundary as possible by varying either tempera- 
ture or interaction strength. Indeed it appears advis- 
able to conduct experiments in systems with the small- 
est order parameter strength for which it is possible to 
reliably visualize the vortex lattice. Both the relatively 
large polarizations and strong interactions required for 
the appearance of j ^ solutions, and the ability to tune 
parameters over wide ranges in atomic systems, demon- 
strate the exceptional potential of tunable cold atom sys- 
tems in the hunt for FFLO vortex latices. The greatest 
obstacle to realization of the FFLO state is likely com- 
petition with phase separated states. We return to this 
point again later. 

The results reported in Fig. (a) can be understood 
qualitatively using quite simple considerations. When 
the temperature is low, weak pairing is expected to be 
dominated by states at the Fermi energy. For that reason, 
the zero-field pairing wavevector on the phase boundary 
is expected to be close to kp-\ — k F [ when T c — > 0, i.e. 
when the interactions are just strong enough to cause 
pairing. Using a small p approximation it follows that 
the pairing wavevector for T c — > is given approximately 
by 



Q 



2k F0 p 



(19) 




The Landau level index at finite fields can be estimated 
by quantizing the pairing wavevector. This gives 



- 1 / k pn a 

FO sc 

(c) Rapid rotation 



FIG. 6: (Color online) Critical T c vs — l/fe^o^sc- (a) 
hClc/epo = 0.02. The curves are for different polariza- 
tions 0.0, 0.1, 0.2, 0.3, 0.4, 0.5 from top to bottom, (b) 
hCl c /EF0 = 0.17. Polarizations are from 0.0 to 0.6. (c) 
HQc/efo = 0.50. Polarizations are from 0.0 to 0.7. 



experimentally when a j ^ vortex lattice occurs. We 
caution, however, that as the temperature drops below 
the critical temperature, different values of j will mix in 
the condensate [3]], |4jj > the j = component will grow 
in weight even if it doesn't have the maximum T c . We 
speculate that the phase transition between finite mo- 
mentum FFLO states and zero-momentum BCS states, 
which occurs at zero field, is replaced in a field by a 
smooth crossover between open and close-packed hexag- 
onal lattices. The best place to search experimentally for 
an exotic vortex lattice is close to the superfluid/normal 



J 



(h 2 Q 2 )/4m e F0 2p 2 



nn c 9 



(20) 



It is easy to check that this equation is quite consistent 
with the numerical results we have obtained. For smaller 
values of hQ c we therefore are confident that even larger 
values of j should occur, although exotic vortex lattice 
may again be confined even more strongly to the region 
close to the phase boundary. For a given value of polar- 
ization, the value of j decreases with increasing interac- 
tion strength because T c moves to higher temperatures, 
reemphasizing the importance of pairing precisely at the 
Fermi energy. 

Fig- ID (b) and (c) show results for systems with larger 
values of h£l c than have been reported in experiments to 
date. One observation is that non-zero j states are less 
likely to occur at large hQ c and appear only at very high 
polarization and strong interactions. This property is 
explained by Eq. lj2U|) . Indeed one can check that the ap- 
pearance of non-zero j values is again consistent with this 
estimate. Other new features that emerge in these fig- 
ures are due mainly to large LL quantization effects. At 
very high rotation frequency [Fig. 0(c)], only the j = 
COM LL is realized. Note that at high temperature, all 
the graphs look similar. (T c decreases monotonically as 
the polarization increases and as the interaction strength 
decreases.) T c is more weakly dependent on the rotation 
frequency. This is because the thermal energy is com- 
parable to or larger than the energy quantization due 
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to rotation. On the other hand, at low temperatures, 
the LL quantization effects become important because 
the particles have one-dimensional densities-of-states for 
each Landau level leading to peaks in pairing (at least 
in this mean-field-theory calculation) when any Landau 
level is just slightly occupied. The non-monotonic density 
of states becomes important when the LL spacing is much 
bigger than the temperature. In this case, we expect 
non-monotonic behavior that is sensitive to the density 
of both hyperfine species; we expect non-monotonic de- 
pendence on polarization and the occasional appearance 
of strong condensates at very large polarizations. Some 
of this non-monotonic behavior is evident in Fig. [B] (c) . 

In Fig. we show the phase boundaries vs. polariza- 
tion and temperature for a series of interaction strengths. 
For slow rotation [Fig. (a)] it is similar to the usual 
BCS-FFLO phase diagram (compare with Fig. |3J). At 
higher rotation frequencies, shown in Fig. 0(b), FFLO 
states are less likely to occur. The transition tempera- 
ture still decreases monotonically as the polarization in- 
creases and above some critical polarization, the normal 
state prevails. At very high rotation frequencies, shown 
in Fig. 0(c), the LL quantization effects become more im- 
portant and we observe reemergence of condensed states 
at around p = 0.4. The difference of the Fermi energies 
at this polarization is exactly equal to the LL spacing and 
the dominant pairing occurs between individual particles 
whose Landau level indices differ by one. 

In Fig.|H|we show an enlargement of the phase diagram 
for the FFLO state, showing also the critical tempera- 
tures for a number of different COM LL index channels 
j in addition to the one with the largest T c . When the 
polarization is small, j = pairing leads to the high- 
est T c ; that is j = is the optimal pairing channel for 
condensation which we denote as j c . As the polariza- 
tion increases, T-! for nonzero j is larger than T^ =0 and 
j c increases with the polarization. This is analogous to 
having an increasing pairing COM momentum with in- 
creasing polarization field in the zero-field case. For a 
given value of Kl Cl non-zero values of j c are more likely 
when interactions are stronger, because the superfluid 
has to be able to withstand the ill effects of polarization 
out to a sufficiently large value of p. If the interaction is 
too weak, no non-zero j pairing can occur and j c is zero. 

In Fig. E] we plot the phase diagram vs. polariza- 
tion and effective interaction space for slow, intermediate, 
and rapid rotations. The critical polarization decreases 
as the interaction strength decreases for weak rotations 
[Fig. El (a)], as seen in experiment The regions la- 
belled FFLO in this figure have j ^ condensates at 
the normal superfluid boundary. Quite generally this be- 
havior occurs only in a small region along the boundary 
between the superfluid and normal state in the regime 
of large polarization and strong interactions. Faster ro- 
tation generally suppresses FFLO states, as emphasized 
earlier, but the superfluid phase can be realized at high 
polarization and weak interaction by tuning the system 
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FIG. 7: (Color online) Polarization vs T c /efo- (a) %Q. c /efo = 
0.02. Curves are for different values of — 1/fcFoOsc- (b) is for 
HQc/efo = 0.17 and (c) is for HQc/efo = 0.50. Dashed lines 
in (a) and (b) shows the T c curves for j = and all the curves 
in (c) corresponds to j — 0. 



such that the Fermi energy mismatch between majority 
and minority species is an integer times the LL spacing, 
and the Fermi energies are close to a quantized LL energy. 
In Fig. El (c) , we see that big peaks occur if these con- 
ditions are met. At zero polarization sf] = £f[ = 1-96, 
in units of LL spacing, and the lowest LL is at 0.5. For 
p = 0.41, £p+ = 2.52 and epi — 1.52 so that the Fermi en- 
ergy difference is exactly the LL spacing and each Fermi 
energy is very close to the LLs. For p = 0.72, e_p| = 2.80 
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FIG. 8: (Color online) Enlarged figures of Fig. |7| (a) for 
— l/fcFoffl sc = (a) 0.6 and (b) 0.8 near the phase boundaries 
between FFLO states and normal fluid. The horizontal axis 
shows T c /efo and the vertical axis is polarization. We cal- 
culate T c for different j'a and determine the optimal j that 
gives the highest T c . 



and efi = 0.80. 



IV. DISCUSSION AND CONCLUSIONS 

In summary, we have derived an equation for the su- 
pcrfluid critical temperature in rotating Fermionic cold- 
atom systems, incorporating Landau level quantization 
effects. Using this equation we have calculated the phase 
boundary between the normal and superfluid phase con- 
sidering pairing in different center-of-mass Landau levels. 
We find that states with higher Landau level condensates 
can occur on the boundary between the normal and su- 
perfluid phase regions in a parameter space that can in 
principle be explored systematically by taking advantage 
of Feshbach resonances and of the ability to create arbi- 
trary degrees of hyperfine state polarization in an atom 
cloud. These FFLO vortex lattice states will have dis- 
tinct vortex lattices [32l l33l l3^ | which should aid their 
identification. High polarization and strong interactions 
are required to realize the FFLO state. At high rotation 
frequency, features that originate from rotational quanti- 
zation effects play an important role and we find that for 
certain parameters the superfluid phase persists to high 
polarization. 

The regime where the FFLO state occurs in rotating 
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FIG. 9: (Color online) Polarization vs —l/kFOdsc- We calcu- 
late T c for polarizations from to I with increment 0.01 and 
choose the largest one that has a finite T c . (a) hQ, c /eF0 = 0.02 
(b) h0. c /e F o = 0.17 and (c) hQ, c /e F0 = 0.50. Solid blue 
curves show phase boundary between normal fluid and su- 
perfluid and dashed blue curves in (a) and (b) show phase 
boundary for COM LL j = 0. 



systems seems accessible to experiment, and hence we be- 
lieve that these exotic vortex structures are observable. 
The greatest obstacle to their observation may be compe- 
tition with states in which the atoms phase separate into 
regions with condensation but no polarization and re- 
gions with polarization but no condensation. We believe 
that FFLO physics would almost certainly occur if phase 
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separation could be suppressed. Phase separation does 
not occur for electrons in a superconducting metal, and 
cannot because of the large Coulomb energy price that 
would have to be paid. One possibility for suppressing 
phase separation in atomic systems with attractive inter- 
actions, is to artificially create the necessary weak but 
long range repulsive interactions by electrically inducing 
dipoles |42( in a pancake shaped |43| , but not necessarily 
quasi-two-dimensional trapped atom system. The typical 
dipole-dipole interaction energy is p 2 /R 3 ~ p 2 n ~ a 2 E 2 n 
where p is the dipolc moment induced by the external 
electric field E, R is the average inter-atom distance, n 
is the density of the atoms and a is the polarizability of 
the atom. If this energy is much smaller than the typi- 
cal atom-atom interaction energy £_Fo(fcFo|<2sc|), then the 
physics on short length scales does not change much. On 
the other hand, if the energy cost of the whole system 
due to the long range dipole interaction when the system 
is phase separated is much larger than the condensation 
energy gain, phase separation can be suppressed. Thus, 
p 2 n 2 V > D(0)A 2 ~ Ne F0 e- 1T / kFa \ a ^ where D{0) is the 
density of states at the Fermi level, N is the number of 
atoms and V is the volume of the system. These con- 
ditions lead to a condition for the external electric field 



e -*/k F0 \a sc \ <<; « <<; fcjro | asc | (21) 

which can be easily satisfied for small kpo\ a sc\- FFLO 
states are most likely expected to occur near the criti- 
cal temperature T c while experimentally observed phase- 
separated states are well below T c . It is known that phase 
separation is less likely at higher temperatures so it could 
be possible to observe FFLO states near T c without ex- 
plicitly suppressing phase separation. 

Finally we mention that peculiar additional interesting 
effects occur because of Landau level quantization if the 
rotation frequency is sufficiently large. Very large rota- 



tion frequencies have been achieved in experiments with 
bosonic atoms [ll. We believe, therefore, that there is no 
fundamental obstacle to approaching the rapid-rotation 
limit with Fermions. Although we have used mean-ficld- 
theory here to study this regime, there is every reason to 
expect unanticipated properties to emerge from strong 
quantum fluctuations and correlations. At sufficiently 
rapid rotations, it should be possible to for the first time 
study the fractional quantum Hall effect in fermion sys- 
tems with attractive interactions |45j . 

This work was supported by the National Science 
Foundation under grant DMR-0606489 and by the Welch 
Foundation. 

APPENDIX: TWO-BODY TRANSITION 
MATRIX AND SCATTERING LENGTH IN 
SYSTEMS WITH ORBITAL MAGNETIC FIELD 

In this appendix, we derive the relation between the 
scattering length and the strength of the delta- function 
like particle-particle interaction in a system with orbital 
magnetic field. The two-body transition operator for 
scattering at energy z is defined by 

f 2B (z) = V + V — —& + ■■■ 
z-H 

= V + V — K^f 2B {z) , (A.l) 
z- H 

where V is the particle-particle interaction and H$ is 
the non-interacting part of the two-body hamiltonian. 
The matrix elements of this transition operator satisfy 
the Lippman-Schwinger equation. Noting that scattering 
conserves the COM motion, we calculate the T-matrix el- 
ements in relative motion Hilbert space with COM LL j 
and COM momenta K y and K z . Notice that the relative 
motion T-matrix does not depend on the y-component 
of the total momentum K y in the Landau gauge. In the 
COM and relative coordinate representation we find 



(N + M- j, k y ,k z \ f 2B (j, K z ; z) \N' + M' - j, k' y ,k' z ) 
= (N + M - j, k y , k z \V\N' + M' - j, k' y X) 



E E 

N",M" k''k" 



(N + M- j, k y ,k z \V\N" + M" - j, k'i, \Bf M " \ 2 (N" + M" - j, k>', k»\V\N' + M' - j, k' yl k' z ) 



Z — £N",K z /2+k'J — £ M",K z /2-k'J 



(A.2) 



Using Eq. we have that, 

(N + M- j, k y ,k z \ 7 L: (j. K, : : ) | A ' + M ' - h' n . /,' ) 

-V 



LyL z 



4'N+M-jO e y^r) < t > N'+M'-j 



1 - 



y 

-V 



\B 



N"M" 12 



Att1 2 L 



v 

B L * N»W'k" Z ~ S N",Kz/2+k'> - £M",Kzl2-k'l 



(A.3) 



r 



For a dilute atomic gas, all the relevant energies are small We are therefore allowed to neglect the energy depen- 
compared to h 2 /mr 2 ? where ry is the interaction range. 
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dence of the two-body T-matrix [33. (Note that the 
energy does not depend on k y .) Hence we have that 

(N + M- j, k yi k z \ f 2B (j, K z -z) \N' + M' ~ j, k' y , k' z ) 
« (0, k y , 0| T 2B (j = 0,K z = 0;z = 0) |0, k' y , 0) 

— M) lr /i_ ;2\ at 11J i2\ 



LyL 



■V 



E 



\Bf M "\* 



4tt1 2 r L z ^> en" k' 
B N",M",k'J 



SM",-k': 



(A.4) 



To extract an expression for the scattering length we put 
the above matrix element equal to the matrix element 
(N + M-j,k y ,k z \ Vp P \N' + M'-j,k' y ,k' z ) of the pseudo- 
potential Vpp(r) = 4:TTa sc h 2 S(r)/m. From this we find 
that 



Anh 2 a s 



1 



1 



4tt1 2 L 



E 



\B, 



NM 1 2 



N~JM~,k z £N ^ + EM >- k * 



(A.5) 
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